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Abstract. We present the CTEQ6HQ parton distribution set which is determined in the general 
variable flavor number scheme which incorporates heavy flavor mass effects; hence, this set pro- 

■ vides advantages for precision observables which are sensitive to charm and bottom quark masses. 
| We describe the analysis procedure, examine the predominant features of the new distributions, 

and compare with previous distributions. We also examine the uncertainties of the strange quark 
distribution and how the the recent NuTeV dimuon data constrains this quantity. 

> 

. Parton distributions functions (PDFs) provide the essential link between the theo- 

retically calculated partonic cross-sections, and the experimentally measured physical 
cross-sections involving hadrons and mesons. This link is crucial if we are to make inci- 



O 

in 



sive tests of the standard model, and search for subtle deviations which might signal new 



physics. The choice of the renormalization scheme is an important issue to address if we 
are to make the most efficient use of our fixed-order perturbation expansion. Separately, 
it is important to know the uncertainty range of the PDFs, and properly fold these into 
O ! the overall uncertainty estimates. We will address both of these issues in turn. 

> 

>< : Global Analysis and the CTEQ6HQ PDFs: 

S' 

The CTEQ6HQ (or C6HQ for short) PDFs[l] are obtained by performing a global 
analysis using the generalized (non-zero quark-mass) MS perturbative QCD framework 
of Refs. [2, 3], which we label the general-mass variable- flavor-number scheme (GM- 
VFNS). When matched to the corresponding hard-scattering cross-sections calculated 
in the same scheme, the combination should provide a more accurate description of the 
precision DIS structure function data, as well as other processes which are sensitive to 
charm and bottom mass effects. 

The C6HQ global fitting follows the same procedure as that of the earlier CTEQ6 
analysis. [4] The data sets used before are supplemented by the HI and ZEUS data sets 
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for the structure function with tagged charm particles in the final state. The F£ data 
sets are quite relevant for this analysis since F 2 C is sensitive to the charm and gluon 
distributions, which are tightly coupled in the generalized MS formalism. 

The C6HQ set is the best fit obtained with these inputs. A broad measure of the 
quality of this fit is provided by the overall % 2 of 2008 for a total number of 1925 
data points (# 2 /DOF = 1.04). This is to be compared to a % 2 of 2037 for 1925 points 
(^ 2 /DOF = 1 .07) in the case of CTEQ6M (or C6M for short). The new C6HQ fit reduces 
the overall j} by 29 out of ~2000 as compared to the C6M fit. The improvement 
of this generalized MS result over the zero-mass MS result is encouraging, since the 
generalized MS formalism represents a more accurate formulation of PQCD. However, 
a difference of % 2 of 29 is within the current estimated range of uncertainty of PDF 
analysis. [4] Therefore, the significance of this difference is arguable. We also note that 
the improvement in j} is spread over most of the data sets: there is no "smoking gun" 
for the overall difference. 

Since perturbative calculations are renormalization scheme dependent, it is important 
to use properly matched hard- scattering cross-sections and PDFs when evaluating fac- 
torized cross-sections for physical applications. This issue is particularly relevant for 
applications involving heavy quarks, since the heavy quark introduces a new mass scale 
which leads to complications of the PQCD formalism. To illustrate this point, we com- 
pare the above results with two possible uses of the PDFs that represent a mis-use of 
PQCD in principle, but occur frequently in the literature in practice, perhaps out of ne- 
cessity. These involve using PDFs obtained in the general-mass scheme convoluted with 
hard- scattering cross-sections (Wilson coefficients) defined in the zero-mass scheme, 
and vice versa. For example, if we use the C6M PDFs which are derived in the ZM-VFN 
scheme with the GM-VFN hard-scattering cross-sections, we obtain a total % 2 of 243 1 
for 1925 data points (# 2 /DOF = 1.26). Conversely, if we use the C6HQ PDFs which are 
derived in the GM-VFN scheme with the ZM-VFN hard-scattering cross-sections, we 
obtain a total % 2 of 2496 for 1925 data points (^ 2 /DOF = 1.30). For the same data sets, 
these mis-matched schemes result in a overall % 2 difference of 420~490 compared to 
the C6HQ set (x 2 of 2008). These are quite large differences relative to the tolerances 
discussed in Refs. [5, 6, 7, 4], and result in obvious discrepancies with some of the pre- 
cision DIS data sets. Clearly, for quantitative applications it is imperative to maintain 
consistency between the PDFs and the hard-scattering cross-sections. 



Comparison with related PDFS: 

The C6HQ and C6M fits provide comparable descriptions of the global QCD data in 
two different schemes. Some of the differences in the PDFs arise purely from the choice 
of scheme. We are particularly interested in the differences for the gluon distribution, 
which will strongly influence the closely correlated charm distribution (since it is gener- 
ated via the g — > cc process). It is also interesting to compare the differences between the 
earlier CTEQ5HQ (C5HQ) distributions with the new C6HQ distributions; differences 
between these PDFs are attributable both to new data, and to minor differences in the 
way the theoretical inputs are implemented. 
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Figure 1. Comparison of the a) gluon and b) strange-quark distributions at Qq = m c = 1.3 GeV for the 
CTEQ5HQ, CTEQ6M, and CTEQ6HQ sets. The axes are scaled to highlight the valence components of 
these distributions. 
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Figure 2. a) Ratio of the 40 CTEQ6M PDF sets compared to the central set for the strange-quark vs x. 
b) Same as previous, with additional sets included (with a more general s(x) parameterization). In both 
figures, guide-lines are drawn at ±5%. 



Fig. la) shows the comparison of the gluon distribution at Qq. Here the difference 
between C5HQ and the CTEQ6 generation of gluon distributions is pronounced. The 
change in this least-well-determined parton distribution is due to the recent precision 
DIS data (most influential in the small x region) in conjunction with the greatly improved 
inclusive jet data from the Tevatron (critical for the medium to large x regions). The 
differences between C6HQ and C6M gluons at large x are due to a combination of 
scheme-dependence, and the inherent uncertainty range of the current analysis. 

Fig. lb) shows the comparison of the strange distributions at the same Qq. The 
noticeable difference between the C5HQ curve and the others, in this case, is largely 
the result of different theoretical inputs: namely, the K parameter which determines the 
ratio of strange to non-strange sea quarks at the initial scale Qq. This K" factor, known 
only approximately, was chosen to be 1/2 both in the CTEQ5 and CTEQ6 analyses, but 
for slightly different values of Qq: 1.0 GeV for C5HQ, and 1.3 GeV for the CTEQ6 
sets. We now look at the uncertainty of the strange quark in further detail; this has 
received increased attention recently since this is an important ingredient in the NuTeV 
determination of singly. 
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Strange-quark PDF and uncertainties: 



Using the set of 40 C6M PDF sets, we can produce a band of of distributions which, in 
principle, should characterize the uncertainty of the s-quark, cf., Fig. 2a). Based on this 
plot, one would be tempted to conclude that the uncertainty on the strange quark is better 
than 5% over much of the x range. Because the parameterization of s(x) in the global 
fit is constrained to be k[u(x) +d(x)] at Qq, the figure actually reflects the uncertainty 
not on s(x) but instead on [ii(x) +d(x)]. This constraint is imposed because none of the 
data in the global analysis directly measures the s(x) distribution; hence, Fig. 2a) is by 
no means a representation of the true uncertainty. 

Recent measurements of the charged-current charm production (vs — > c/i — > /i ± /i =F X) 
from CCFR and NuTeV [8, 9] have the potential to determine {s(x) 1 s(x)} with much 
more precision. In Fig. 2b) we again show the 40 C6M PDF sets combined with a 
number of additional fits that relax the parameterization of s(x) . This is not an exhaustive 
collection designed to span the full parameter space, but rather an illustration that the 
implied uncertainty of Fig. 2a) is much too conservative; the dimuon data will go a long 
way toward improving our knowledge of the strange PDF. An accurate determination 
of the strangeness of the proton, as well as the strangeness asymmetry [s(x) —s(x)] will 
have important implications for a number of measurements, including the NuTeV sin 0w 
measurement. 



Concluding Remarks: 

The C6HQ PDFs presented here complement the CTEQ6 sets by providing distri- 
butions which can be used in the generalized MS scheme with massive partons. This 
analysis includes the complete set of NLO processes including the real and virtual quark- 
initiated terms. 

While the zero-mass parton scheme is sufficient for many purposes, the fully massive 
scheme can be important when physical quantities are sufficiently sensitive to heavy 
quark contributions. This is evident when comparing the C6HQ and C6M fits to the 
mis-matched sets where the precise DIS data from HERA highlights the discrepancies. 

The C6HQ fits also provide the basis for a series of further studies involving more 
quantitative analysis of strange, charm, and bottom quark distributions inside the nu- 
cleon. For example, the C6HQ PDFs are necessary for a consistent analysis of resumed 
differential distributions for heavy quark production. [10] Using the full range of data 
from both the charged and neutral current processes, these distributions can reduce the 
uncertainties in the calculations; hence, they have significant implications for charm 
and bottom production, and can help resolve questions about intrinsic heavy quark con- 
stituents inside the proton, the AxFt, structure function, and the extraction of sin 6yy- 



Heavy Quark Parton Distribution Functions 



2nd February 2008 



4 



Acknowledgments: 



F.I.O. acknowledges the hospitality of Fermilab and BNL, where a portion of this 
work was performed. This work was supported by RIKEN, BNL, the U.S. DoE under 
grant DE-AC02-98CH10886 & DE-FG03-95ER40908, and the Lightner-Sams Founda- 
tion. 



REFERENCES 

1. S. Kretzer, H. L. Lai, F. I. Olness, and W. K. Tung, Phys. Rev., D69, 114005 (2004), hep-ph/ 
0307022. 

2. M. A. G. Aivazis, J. C. Collins, F. I. Olness, and W.-K. Tung, Phys. Rev., D50, 3102-31 18 (1994), 

hep-ph/9312319. 

3. W.-K. Tung, S. Kretzer, and C. Schmidt, J. Phys., G28, 983-996 (2002), hep-ph/0110247. 

4. J. Pumplin, et al., JHEP, 07, 012 (2002), hep-ph/02 01195. 

5. A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thome, J. Phys., G26, 663-665 (2000). 

6. A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thome, Eur. Phys. J., C28, 455-473 (2003), 

hep-ph/0211080. 

7. A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thome, Eur. Phys. J., C35, 325-348 (2004), 

hep-ph/0308087. 

8. A. O. Bazarko, et al., Z. Phys., C65, 189-198 (1995), hep-ex/ 94 60 07. 

9. M. Goncharov, et al., Phys. Rev., D64, 1 12006 (2001), hep-ex/ 102 04 9. 

10. P. M. Nadolsky, N. Kidonakis, F. I. Olness, and C. P. Yuan, Phys. Rev., D67, 074015 (2003), 

hep-ph/0210082. 



Heavy Quark Parton Distribution Functions 



2nd February 2008 



5 



